Strange Particle Production 
from Quark Matter Droplets 



Klaus WERNER ^ and Michael HLADIK 



Institut fur Theoretische Physik, Universitdt Heidelberg, Germany 



Abstract. We recently introduced new methods to study 
ultrarelativistic nuclear scattering by providing a link between 
the string model approach and a thermal description. The string 
model is used to provide information about fluctuations in energy 
1— ' density. Regions of high energy density are considered to be quark 

,-0 matter droplets and treated macroscopically. At SPS energies, we 

ry^ find mainly medium size droplets — with energies up to few tens 

of GeV. A key issue is the microcanonical treatment of individual 
quark matter droplets. Each droplet hadronizes instantaneously 
according to the available n-body phase space. Due to the huge 
O^ number of possible hadron configurations, special Monte Carlo 

techniques have been developed to calculate this disintegration. 
We present results concerning the production of strange particles 
from such a hadronization as compared to string decay. 
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I. INTRODUCTION 



Studying nuclear collisions at ultrarelativistic energies (E'cms/nucleon ^> 
1 GeV) is motivated mainly by the expectation that a thermalized system 
of quarks and gluons (quark-gluon plasma) is created [1]. There are essen- 
tially two directions for modelling such interactions: dynamical and thermal 
approaches. The former ones refer to string models [2-7] or related methods 
[8], supplemented by semihard interactions at very high energies [9-12]. Here, 
a well established treatment of hadron-hadron scattering, based on Pomerons 
and AGK rules [13], is extended to nuclear interactions. Thermal methods 
[14-19] amount to assuming thermalization after some initial time r , with 
evolution and hadronization being mostly based on ideal gas assumptions. 
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FIG. 1. High density regions in the x-£ plane for a typical event. 

Both methods have serious theoretical drawbacks. Even for nuclei as 
light as sulfur the string models produce particle densities that high that 
the hadrons are overlapping. So the independent string model is certainly 
too simplistic, and also considering secondary interactions as binary collisions 
among hadrons can theoretically not be justified. On the other hand it is also 
not realistic to consider a homogeneous plasma occupying the whole available 
volume, as assumed in thermodynamic models. To illustrate this, we show 
in fig. 1 a typical event of a string model simulation. We consider a central 
S+S collision at 200 GeV (E cms /A 10 GeV), the transverse coordinates 
being x and y, the longitudinal one (= beam direction) being z\ it is useful 
to consider space-time rapidity ( = 0.51n(t + z)/(t — z) rather than z. In the 
figure, the hatched regions represent high energy density (e > e c = 1 GeV) in 
the x — ( plane (y = 0). We find a couple of intermediate size regions of high 
energy density, representing rest masses of few GeV up to few tens of GeV. 
This demonstrates that neither the independent string model is correct, since 
these high density regions cannot possibly be treated as independent hadrons 
nor the thermal approaches, since we do not have one big high density object 
but rather a couple of medium size objects in addition to plenty of ordinary 
hadrons and resonances in particular in the periphery. 

In this paper we sketch our new approach [20,21], which is more realistic 
than the string model and more realistic than thermal treatments, providing a 
link between the two. We then show first results, where we compare rapidity 
and m t spectra from QM droplet hadronization with the corresponding dis- 
tributions from string decay. Results from full nucleus-nucleus collisions will 
be published soon. 
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II. THE NEW APPROACH 



In this section, we review the basic features of our approach, which is 
meant to be a link between the "bare" string model and thermal treatments. 

Based on the string model, we first determine connected regions of high 
energy density (e > e c , for a given e c ). These regions are referred to as 
quark matter (QM) droplets. For such regions, the initially produced hadrons 
serve only as a mean to produce the proper fluctuations in the energy density. 
Presently, a purely longitudinal expansion of the QM droplets is assumed. 
Once the energy density falls beyond e c , the droplet D decays instantaneously 
into an n-hadron configuration K = {hih 2 ■ ■ -h n } with a probability propor- 
tional to 0, with representing the microcanonical partition function of an 
n-hadron system. Due to the huge configuration space, sophisticated methods 
of statistical physics [22,23] have to be employed to solve the problem without 
further approximations. 

Our hadronization scenario is referred to as "microcanonical hadron gas 
(MHG) scenario". It is certainly not the only one and probably not the most 
realistic one. However, we start with this scenario for a couple of reasons: the 
hadron gas scenario is a benchmark, widely used in the literature (in a sim- 
plified fashion though); the MHG scenario can be solved exactly; for massless 
hadrons even an analytical treatment is possible, providing very useful tests 
for the complicated numerical procedures. After having gained experience in 
the techniques to solve the MHG scenario, we plan to investigate alternatives 
as well. So the purpose of this paper is not so much to promote this particular 
scenario, but rather to show how a dynamical and a statistical treatment can 
be combined. 

The first stage of our approach is the identification of high energy density 
regions, based on the string model, which is already discussed elsewhere [20]. 
Due to the empirically found correlation, y = £, between the average rapidity y 
and space-time rapidity £, a hypersurface 7i T of constant proper time r may be 
introduced, in the central region simply defined by t 2 -z 2 = t\ Appropriate 
coordinates on 7i T are the space-time rapidity ( = 0.5 ln(t + z)/(t — z) and the 
transverse coordinates x and y. After having used the string model (VENUS 
5.08) to get complete information of hadron trajectories in space and time, 
we may now, for given r, determine energy densities on 7i T and thus locate 
high density regions on 7i T (with e > e c ), as shown in figure 1 for a typical 
example. 

High density regions are considered as QM droplets, presently it is assumed 
that they expand purely longitudinally. Whenever other clusters or hadrons 
cross their way, the two objects fuse to form a new, more energetic cluster. 
Due to the expansion, the energy density of a cluster will at some stage drop 
below e c , which causes an instantaneous decay. 
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We employ the "micro-canonical hadron gas (MHG) scenario" for the had- 
ronization: the probability of a droplet D — characterized by the invariant 
mass E, the volume V, and the flavour Q = (Q u , Q d , . . .) — to decay into a 
hadron configuration K = {hi, . . . , h n } of hadrons hi is given as 

probp -> K) ~ tt(K) , (f) 

with 0(A") being the microcanonical partition function of an ideal, relativistic 
gas of the n hadrons hi, 



0,(K) — C V ol Cdeg C'ident 4> 



Wl 
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^vol = (C) s 3n , Cdeg = II 9i , Cident = II ~ f - ( 3 ) 

) 1=1 aes 

Here, Cd eg accounts for degeneracies ((/; is the degeneracy of particle i), and 
Cident accounts for the occurence of identical particles in K (n a is the number 
of particles of species a). The last factor 

~ n 

4> = j n <?Pi s ( e - Ee o s &Pi) 8 Qt x qi (4) 



8 = 1 



is the so-called phase space integral, with e 8 - = yrrq + pf being the energy 
and f>i the 3-momentum of particle i. The vector qi = (g", qf, . . .) represents 
the flavour content of hadron i. The expression eq. (4) is valid for the centre- 
of-mass frame of the droplet D. 

We have to define a set S of hadron species; we include the pseudoscalar 
and vector mesons (it, K,rj,r]', p, K* ,cu,c/>) and the lowest spin-| and spin- 
| baryons (N, A, E, S, A, E*, £*, 0) and the corresponding antibaryons. A 
configuration is then an arbitrary set {hi, . . . , h n } with hi £ S . 

We are interested in droplet masses from few GeV up to 10 3 GeV, corres- 
ponding to particle numbers n = \K\ between 2 and 10 3 . So we have to deal 
with a huge configuration space, which requires to employ Monte Carlo tech- 
niques, well known in statistical physics. The method at hand is to construct 
a Markov process, specified by an initial configuration K , and a transition 
probability matrix p(K t — > Kt+i)- In generating a sequence K , Ki, K 2 , ■ ■ ., 
two fundamental issues have to be payed attention at: 

• initial transient: starting usually off equilibrium, it takes a number of 
iterations, 7 eq , before one reaches equilibrium; 

• autocorrelation in equilibrium: even in equilibrium, subsequent config- 
urations, K a and K a+ i, are correlated for some range I au to of i. 
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In general, both 7 eq and I au to should be as small as possible. 

We are going to proceed as follows: for a given droplet D with mass A, 
volume V, and flavour Q, we start from some initial configuration K 0} and 
generate a sequence K 0} A"i, . . . , Ki , with 7 eq being sufficiently large to have 
reached equilibrium. If we repeat this procedure many times, getting con- 
figurations Kj^,Kj^,..., these configurations are distributed as fi(A"). So 
for our problem, we have only to deal with the initial transient, not with the 
autocorrelation in equilibrium. We have to find a transition probability p such 
that it leads to an equilibrium distribution fi(A"), with the initial transient 
A. q being as small as possible. 

Sufficient for the appropriate convergence to fi(A") is the detailed balance 
condition, 

n(K a ) p(K a ^ K h ) = n(K b ) p(K b ^ K a ) , (5) 

and ergodicity, which means that for any K a , Kb there must exist some r with 
the probability to get in r steps from K a to Kb being nonzero. Henceforth, 
we use the abbreviations 

fi a := VL(K a ); p ab := p(K a -> K h ). (6) 

Following Metropolis [22], we make the ansatz 

Pab = Wab U a b , (7) 

with a so-called proposal matrix w and an acceptance matrix u. Detailed 
balance now reads 



U a b _ Qb W ba 
U ha a W a b 

which is obviously fulfilled for 



f8) 



Uab = F[- , 9 

\\la WabJ 

with some function F fulfilling F(z) / F(z~ 1 ) = z. Following Metropolis [22], 
we take 

F(z) = min(2,l) . (10) 

The power of the method is due to the fact that an arbitrary w may be chosen, 
in connection with u being given by eq. (9). So our task is twofold: we have 
to develop an efficient algorithm to calculate, for given K, the weight 0(A"), 
and we have to find an appropriate proposal matrix w which leads to fast 
convergence (small 7 eq ). The first task can be solved, a detailed description 
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will be published soon [24]. In the following we discuss about constructing an 
appropriate matrix w. 

Most natural, though not necessary, is to consider symmetric proposal 
matrices, w a b = Wf, a} which simplifies the acceptance matrix to u a b = 
F(ftb/rt a ). This is usually referred to as Metropolis algorithm. Whereas for 
spin systems, it is obvious how to define a symmetric matrix w, this is not so 
clear in our case. We may take spin systems as guidance. A configuration K 
is per def. a set of hadrons {hi, . . . , h n } with the ordering not being relevant, 
so {7r°,7r°,p} is the same as {p, 7r°,7r }. We introduce "microconfigurations" 
to be sequences {hi, . . . , h n } of hadrons, where the ordering does matter. So 
for a given configuration K a = {hi, . . . ,h n } there exist several microconfig- 
urations K a j = {/^(i), . . . , /i7rj(n)} 5 with 7Tj representing a permutation. The 
weight of a microconfiguration is 




with n a being the number of hadrons of type a. Taking for example 
K = {p, 7T°,7r }, there are three microconfigurations {p, 7r°,7r }, {ir°,p, ir } 
and {7r°,7r°,p}, with weight 0(A")/3. 

So far we deal with sequences {hi, . . . , h n } of arbitrary length n, to be 
compared with spin systems with fixed lattice size. We therefore introduce 
"zeros", i.e. we supplement the sequences {hi, . . . , h n } by adding L — n zeros, 
as {hi, . . . , h n , 0, . . . , 0}, to obtain sequences of fixed length L. The zeros may 
be inserted at any place, not necessarily at the end. Therefore the weight of 
a microconfiguration K a j with zeros relative to the one without, K a j, is one 
divided by the number of possibilities to insert L — n zeros, so from eq. (11) 
we get 

= ^ (n ™A nl[L ~ n)l n(K a ) . (12) 

We now have the analogy with a spin system: we have a one-dimensional 
lattice of fixed size L, with each lattice site containing either a hadron or 
a zero. Henceforth, we use for microconfigurations with zeros the notation 
K aj = {hi, . . . , hzj} with hi being a hadron or zero. 

Since from now on we only consider microconfigurations with zeros (K a j) 
rather than configurations (K a ), we are going to write K a instead of K a j, 
keeping in mind that a represents a double index, and say "configuration" 
rather than "microconfiguration with zeros". The advantage is that we can 
use the above formulas specifying the Metropolis algorithm without changes. 

We are now in a position to define a symmetric proposal matrix w(K a — > 
K b ), with K a = {hi, ...,h a L } and K b = {h\, . . . , h b L }, as 



6 




1000 



1500 2000 
iterations 



FIG. 2. Multiplicity versus the number of iterations. 
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with 



WW **W«»W»0 . (14) 

where V{h a i h a j j ) is the set of all pairs (hihj) with the same total flavour as the 
pair (hfh'j). The symbol \V\ refers to the number of pairs of V. The term {} 
in eq. (13) makes sure that up to one pair all hadrons in K a and Kb are the 
same, the term 2/ L(L — 1) is the probability to randomly choose some pair of 
lattice indices i and j. So our proposal matrix amounts to randomly choosing 
a pair in K a , and replacing this pair by some pair with the same flavour, 
with all possible replacements having the same weight. The proposal matrix 
is obviously symmetric, since v is symmetric (the symmetry of v is crucial!). 
We have now fully defined an algorithm, which due to general theorems will 
converge, but how fast, i.e., how large is 7 eq ? Considering particle ratios, 
like n v o /n jr + , we find immediately that we have a very slow convergence, 
so igq is too large for the method to be of practical importance. This is 
obvious, since the method is not very democratic: flavourless particles like 
7T°, p° or also zeros are much more frequently proposed than all the rest. This 
shortcoming can be fixed by defining w such that two pairs are exchanged 
rather than one, the first pair being replaced by a completely arbitrary pair, 
the second one by some pair to guarantee flavour conservation. In addition it 
is necessary to weight the "zeros" differently than the hadrons. This improved 
method violates the symmetry of w, however, the asymmetry w a b/wb a can be 
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FIG. 3. Comparing hadronization of quark matter droplets (solid) and 
string decay (dotted): Rapidity spectra (left) and transverse mass distribu- 
tions (right) of protons (top) and antiprotons (bottom). 

calculated and properly taken into account. Further details of the "asymmetric 
algorithm" will be published elsewhere [24]. 

The asymmetric method converges quite fast. As a check, we consider 
massless hadrons, where analytical results exist. In fig. 2 we plot the iterated 
total multiplicity N for a droplet with E = 10 GeV and V = 10 fm 3 , compared 
to the average multiplicity N from the analytical calculation (dashed line). 
The "equilibration time" 7 eq is roughly given as 

Zeq/iVwlO, (15) 

with the initial configuration being two 7r 's. 

III. RESULTS 

In the following we compare hadron production from QM droplets — as 
discussed in the previous section — with string decay. We consider droplets 
with E = 10 GeV, V = 10 fm 3 and uuo?-flavour and correspondingly u-ud 
strings with an energy of 10 GeV. We are going to present rapidity (y) and 
transverse mass (m t ) spectra for different hadrons, in particular strange ones. 
In fig. 3, we show results for protons and antiprotons. Here — as in the 
following figures — the solid histograms represent QM droplets, the dashed 



8 



10 r 



; 1 GeV - 


> A 




E 


-1 


F 1 GeV - 


> A 






VENUS 5.08 


s 


10 








VENUS 5.08 


- 




















J \ 0M 




e" 










QM - 






g ... 




10 








, string ... 










-3 
10 




















-4 

10 












nl, , i , 












/ = o 






\ , , , i 














I 







4 -2 





2 

rapidity y 


- 1 GeV — > 


A 








VENUS 5.08 






QM - 






string ... 




, , I 


1 



1 1.2 1.4 1.6 1.8 2 
tronsv. mass m, 



: 1 GeV -> A 








VENUS 5.08 






QM - 






string ... 




y - : : 





2 4 
rapidity y 



1 1.2 1.4 1.6 1.8 2 
transv. mass m. 



FIG. 4. Comparing hadronization of quark matter droplets (solid) and 
string decay (dotted): Rapidity spectra (left) and transverse mass distribu- 
tions (right) of lambdas (top) and antilambdas (bottom). 

ones strings. For protons from strings, we clearly observe a "leading particle 
effect", i.e. a maximum of the y distribution at forward rapidities. This is 
due to the fact that the fast forward ud diquark is constituent of the observed 
proton. In case of antiprotons, we have a "real" production process with no 
string (or droplet) constituents appearing in the hadron. We clearly observe 
more antiprotons from QM droplets compared to string decay, the m t spectra 
in the QM case are somewhat flatter. The (anti)lambda spectra, shown in 
fig. 4, are quite similar to the (anti)proton case, but the A enhancement 
(of QM compared to string decay) is large, roughly a factor of ten. For the 
(anti)cascades, presented in fig. 5, we find even more enhancement, about a 
factor of fifty. 

The baryon-antibaryon (BB) production may be summarized as follows: 
whereas BB pairs are rarely produced from string decay, QM droplets provide 
considerably larger rates. The enhancement increases with increasing strange- 
ness content (the string results drop dramatically). Such large enhancement 
(of 10 or 50) will not show up in nucleus-nucleus collisions, since here we have 
a mixture of hadrons from droplets and those from string decay (formed in 
dilute areas). 

We now turn to mesons. The pion rapidity spectra from strings, shown 
in fig. 6, are much broader than the spectra from droplet decay, the absolute 
numbers (integrals) being not too different though. Also the m t spectra are 
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FIG. 5. Comparing hadronization of quark matter droplets (solid) and 
string decay (dotted): Rapidity spectra (left) and transverse mass distribu- 
tions (right) of cascades (top) and anticascades (bottom). 

quite similar. For the kaons (K + ) } as shown in fig. 7, we find some enhance- 
ment, much less though as compared with A and S. 

So comparing droplet with string decay, we find a substantial enhancement 
of antibaryons — the more strangeness the better. This is due to the fact that 
BB production is strongly suppressed in string decay, with no such suppression 
acting in droplet decay. 



IV. SUMMARY 

We have presented our new method to model ultrarelativistic scattering, 
where one attempts to link the string model and statistical approaches. The 
string model is used to locate high density regions, which are then referred 
to as quark-matter droplets and treated macroscopically. We presented first 
results which showed that antibaryons, in particular strange ones, are much 
more frequent in droplet decay as compared to string decay. This will lead 
to an enhanced antibaryon rate in nuclear collisions as compared to proton- 
proton. 
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FIG. 6. Comparing hadronization of quark matter droplets (solid) and 
string decay (dotted): Rapidity spectra (left) and transverse mass distribu- 
tions (right) of 7r + 's (top) and 7r~'s (bottom). 
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